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I. INTRODUCTION 

Non equilibrium phenomena are on the agenda in modern statistical physics, soft condensed matter and bio- 
physics. Open systems driven far from equilibrium are ubiquitous. A classical case is driven Navier Stokes turbulence, 
other cases are driven lattice gases, growing interfaces, growing fractals, etc. Unlike equilibrium physics where the 
Boltzmann-Gibbs scheme applies, the ensemble is not known in non equilibrium. Here the problem is defined in terms 
of a numerical algorithm, a master equation, or a Langevin equation. 

An interesting class of non equilibrium systems exhibit scale invariance. One example is diffusion limited aggregation 
(DLA) driven by the accretion of random walkers yielding a growing scale invariant fractal with dimension D « 1.7 in 
2d. Another case is a growing interface driven by random deposition or propagating in a random environment. Here 
the width of the growing front w(L,t) conforms to the dynamical scaling hypothesis w(L,t) = F> f{tj L z ), where L is 
the size of the system and £ and z scaling exponents; Q characterizing the 

roughness and z describing the dynamical crossover to the stationary profile [l[ . 

In the present paper we focus on the Kardar-Parisi- Zhang (KPZ) equation which describes an intrinsic non 
equilibrium problem and plays the same role as the Ginzburg-Landau functional in equilibrium physics. The KPZ 
equation was introduced in 1986 in a seminal paper by Kardar, Parisi and Zhang 0. It has the form 

dh A 

— = v V 2 h+-VhVh-F + i 1 , (w)(f,t) = A5(r)6(t), (1.1) 

and purports to describe non equilibrium aspects of a growing interface; see Refs. [HQ. Here h(r,t) is the height 
of an interface at position r and time t, the linear diffusion term v'S/ 2 h, characterized by the diffusion coefficient u, 
represents a surface tension, the nonlinear growth term (\/2)VhS7h, characterized by A, is required to account for 
the lateral growth, F is an imposed constant drift, and the random aspects, i.e., the random deposition of material or 
the random character of the medium, are encoded in the noise rj(r, t). The noise is assumed to be locally correlated 
in space, and time, its strength characterized by A. 

Despite its simple form the KPZ equation is difficult to analyze and many aspects remain poorly understood [J [H, HI • 
Apart from its intrinsic interest the KPZ equation is also related to fundamental issues in turbulence and disorder. 
Introducing the local slope field u = V/i the KPZ equation takes the form of a Burgers equation driven by conserved 
noise, 

^ - X(u ■ V)u = vV 2 u + V?y. (1.2) 

In the noiseless case setting A = — 1 and regarding u as a velocity field, the non linear term appears as the convective 
term in the Navier-Stokes equation and Eq. (jl.2p has been used to model aspects of turbulence. In the Id case the 
relaxation of the velocity field takes place subject to a transient pattern formation composed of domain walls and 
ramps with superimposed diffusive modes. In the driven case Eq. (|1.2[) was studied earlier in the context of long time 
tails in hydrodynamics On the other hand, applying the non linear Cole-Hopf transformation the KPZ equation 
maps to the Cole-Hopf equation (CH) 

dw o -» ^ n A 2v 

-7— = vSJ u - — wF + — wr], h=—mw, (1.3) 
ot 2v 2v A 
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a linear diffusion equation driven by multiplicative noise. The CH equation has a formal path integral solution 0, Q 
which can be interpreted as an equilibrium system of directed polymers (DP) with line tension 1/4^ in a quenched 
random potential 77; a model system in the theory of disorder which has been studied using replica techniques 0, IE). 

The KPZ equation lives at a critical point and conforms to the dynamical scaling hypothesis. For the height 
correlations we have 

(hh)(r,t) = r 2C F(t/r z ). (1.4) 

Here C, z, and F are the roughness exponent, dynamic exponent, and scaling function, respectively. To extract scaling 
properties the initial analysis of the KPZ equation was based on the dynamic rcnormalization group (DRG) method, 
previously applied to dynamical critical phenomena and noise-driven hydrodynamics [3]. An expansion in powers of 
A in combination with a momentum shell integration yield to leading order in d — 2 the DRG equation dg/dl = /3(g), 
with beta-function B(g) = (2 — d) (7 + const. g 4 . Here g = AA 2 /^ 3 is the effective coupling strength and I the logarithmic 



scale parameter [HQ. The emerging DRG phase diagram is depicted in Fig. [TJ 
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Figure 1: DRG phase diagram for the KPZ equation to leading loop order in d — 2. In d = 1 the DRG flow is towards the 
strong coupling fixed point C, = 1/2,2 = 3/2. Above the lower critical dimension d = 2 there is an unstable kinetic transition 
line, separating a rough phase from a smooth phase. 

Before discussing the phase diagram we note two further properties of the KPZ equation. First, subject to a Galilean 
transformation the equation is invariant provided we add a constant slope to h and adjust the drift F, i.e., 

r^r-\u°t, h^h + v°-r, F -> F + (\/2)(u°) 2 ; (1.5) 

note that the slope field u and the diffusive field w transform like u — > u + iP and w — » w exp[(A/2z/)u° -r\, respectively. 
The Galilean invariance implies the scaling law 

( + z = 2, (1.6) 

relating £ and z Q; the Galilean invariance is a fundamental dynamical symmetry specific to the KPZ equation, 
delimiting the universality class. Second, a fluctuation-dissipation theorem is operational in Id since the stationary 
Fokker-Planck equation admits the explicit solution Q 



Po(h) oc exp 



-(f/A) / dx{Vh) 



(1.7) 
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This distribution shows that the slope u = Vh fluctuations are uncorrected and that the height field h = f x udx' 
performs a random walk. Consequently, from Eq. (jl.4p we infer the roughness exponent C, = 1/2 and from the scaling 
law (|1.6j) the dynamic exponent z = 3/2. In other words, the scaling exponents associated with the strong coupling 
fixed point are exactly known in Id, see Fig. [TJ moreover, results for the scaling function can be obtained by loop 
expansions Q. 

The lower critical dimension is d = 2. Below d = 2 the DRG flow is towards the strong coupling fixed point. Above 
d = 2 the DRG equation yields an unstable kinetic phase transition line as indicated in Fig. [TJ For A below a critical 
coupling strength A c the DRG flow is towards A = 0, corresponding to the linear Edwards- Wilkinson (EW) equation, 
the KPZ equation for A = 0, yielding z — 2 and £ = (2 — d)/2, note that the scaling law is not operational for A = 0; 
for A > A c the DRG flow is towards a non perturbative strong coupling fixed point. 

A DRG analysis to all order in d — 2 yields z = 2 and C = on the transition line and a singularity at the upper 
critical dimension d_un pcr = 4. Mode coupling techniques also give d U ppcr = 4, whereas a directed polymer analysis 
yields d U p por w 2.5 0. 

In order to disentangle the properties of the KPZ equation, both regarding scaling and otherwise, there is a need 
for alternative methods. In the present paper we summarize a non perturbative weak noise approach which we have 
pursued in recent years Q. The working hypothesis here is to focus on the noise strength A as the determining 
parameter rather that the non linearity A and ensuing DRG analysis or mapping to DP. For A = the interface 
relaxes subject to a transient pattern formation; for A w the interfaces initially decays but is eventually driven into 
a stationary fluctuating state; the crossover time diverging for A — > 0. In the next section we summarize the scheme 
that will allow us to access the weak noise regime in a non perturbative fashion. 



II. WEAK NOISE SCHEME 



The weak noise scheme takes as its starting point a generic Langevin equation 

dx 



dt 



= -F(x)+r,, ( W )(t) = A6(t), 



(2.1) 



determining the problem we wish to analyze. Here a; is a multi-dimensional random variable, F(x) a general non 
linear drift, and r\ additive white noise correlated with strength A. In order to implement a weak noise approximation 
we consider the equivalent Fokkcr-Planck equation for the distribution P(x, t) 



BP 1 a2 S 2 P 
A — = -A 

dt 2 dx 2 



(2.2) 



Interpreting Ad/dx as a momentum operator, P as an effective wave function, and A as an effective Planck constant, 
Eq. (|2.2p has the form of an imaginary time Schrodingcr equation and it is natural to introduce the well-known WKB 
or eikonal approximation 



Pocexp[-S/AJ. 



(2.3) 



To leading order in A the action then obeys a principle of least action 5S = as expressed by the Hamilton- Jacobi 
equation dS/dt + H(x,p) — with canonical momentum p = dS/dx. The Hamiltonian (energy) takes the the form 



H 



1 



1 



p 2 - pF(x) = -p[p - 2F(x)} 



yielding the coupled equations of motion 

dx dp 

Tt=- F + P > Tt 

The action associated with an orbit from x\ to x in time T is 



S(x\ — > x, T) 



x,T 



dt 



dx 
'"eft 



H 



dF_ 

dx 

1 

~ 2 



dtp(tf 



(2.4) 



(2.5) 



(2.6) 



£Ci,0 



The weak noise recipe is clear. We solve the equations of motion (|2.5p and identify an orbit from x\ to x in time T 
with p as an 'adjusted' variable. The orbits lie on constant H manifolds. Evaluating the action S for a specific orbit 
the WKB approximation (|2.3[) yields the transition probability from x\ to x in time T. Assuming F(x) oc x for small 
x the phase space has the generic structure depicted in Fig. [21 

The present variationally based weak noise scheme, dates back to Onsager. In more recent formulations it corre- 
sponds to the saddlepoint contribution (optimal path) in the functional Martin-Siggia-Rose scheme, see Refs. [1,0]- 
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Figure 2: Canonical phase space plot. The finite-time orbit from the initial xl to the final x in transition time T lies on 
the energy manifold H 7^ 0. In the long time limit the orbit migrates to the zero-energy manifold composed of a transient 
submanifold for p — 0, corresponding to the noiseless case, and a stationary manifold for p — 2F, corresponding to the noisy 
case. The submanifolds intersect in the saddle point (SP), determining the Markovian behavior. 



III. GROWTH MODES 



The weak noise scheme applies directly to the KPZ equation in the Cole-Hopf formulation (|1.3p . Extending the 
scheme in order to incorporate multiplicative noise, see Refs. the WKB scheme yields the equations of motion, 
action, and distribution 

^ = v[V 2 w - k 2 w] + k 2 w 2 p, ^ = -v[V 2 w-k 2 w]-klpw 2 , (3.1) 

/w,T 
dfdt(wp) 2 , P(w,T) <xexp[-S(w, T) /A], (3.2) 

with parameters k 2 = \F/2v and fco = \/2v. On the transient and stationary manifolds p = and p oc w the 
equations of motion reduce in the static case to the diffusion and non linear Schrodinger equations 

V 2 w = k 2 w, (3.3) 
V 2 w = k 2 w - k^w 3 , (3.4) 

admitting localized spherically symmetric solutions w+(r) oc exp(kr) and W-(r) oc exp(— kr) for large r, respectively. 
In terms of the height field h = k^ 1 \nw and slope field u = V/i we obtain the fundamental localized static growth 
modes for large r 

h±(r) oc ±— r and u±(r) oc ±- . (3-5) 

ko k r 

The growth modes are characterized by the amplitude or charge k oc F 1 ^ 2 , determined by the imposed drift F in the 
KPZ equation. The growth mode with positive charge h oc exp(fcr/fco), u oc r/r lives on the transient 'noiseless' p = 
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Figure 3: We depict the static growth modes in the height and slope fields in 2d. 



manifold and corresponds to a cone or dip in h and a constant positive monopole in u. The mode carries zero action, 
zero energy H = — S/T=0, and zero momentum II = J drw'Vp = 0; the distribution P oc cxp(— S/A) associated 
with the mode is of 0(1). The growth mode with negative charge h oc exp(— kr/ko), u oc — r/r is associated with the 
stationary 'noisy' manifold p oc w and corresponds to an inverted cone or tip in h and a constant negative monopole 
in u. This modes carries a finite action 

S oc T(iVfc ) 2 fc 4-d , ( 3 - 6 ) 

yielding the distribution P oc exp[— S/A]. 

In the Id case the growth modes correspond to right hand and left hand domain walls, u oc ±tanh(for). The right 
hand domain wall is the viscosity broadened shock wave in the deterministic Burgers equation; the left hand domain 
wall is 'noise induced', in the present scheme characterized by a finite p, see Refs. @. In Fig. [3] we have depicted the 
static growth modes in the height and slope fields. 

Applying the Galilean transformation (|1.5p we can boost the static modes and obtain the propagating growth modes 

h±(r,t) = ±—]r + Xu a t\+u -r, u±(r, t) = ± — — — — — + u Q . (3.7) 
fco ho \r + \uot\ 

The moving localized growth modes are the fundamental elementary excitations incorporating the non linear aspects 
of the KPZ equation. In a quantum field theory context the growth modes correspond to instantons or solitons. 



IV. STOCHASTIC PATTERN FORMATION 



By means of the propagating localized growth modes (|3.7p we construct a global solution of the field equations 
J) . The Galilean invariance ()1 .5|) determines the matching of the modes and we obtain the dilute network solution 

i i ' ^ ' ' 

nit) -n(t) 



i(t) = -2uJ2h ^.i J . , , n(t)= / Zi(t')dt> + n(0)- (4-2) 
\n(t)-ri{t)\ J 



Here ki is the assignment of charges and r 2 (0) the initial positions. This network construction correspond to a multi- 
instanton solution in quantum field theory. As time evolves the modes propagate and the velocities adjust to constant 
values given by the self-consistent equation 

^ = -2z,Vfe iT $^! T . (4.3) 

Superimposed on the propagating network of growth modes is a spectrum of linear extended diffusive modes with 
dispersion uj = vk 2 , following from a linear analysis of the field equations (|3.f I) . 

For large r the slope field u ~ (^/^)X^j ^ n order to ensure the boundary condition of a flat interface at large 
distances we impose the neutrality condition £^ fcj = 0; note that this boundary condition still allows for a local 
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Figure 4: We show a dipole configuration composed of two paired monopoles of opposite charges. In a) we depict the dipole 
mode for the height field in 2d; in b) we show the two-domain wall configuration for the slope field in Id. 

offset of h and the propagation of facets or steps. During time evolution the dynamical network propagates across 
the system. Imposing periodic or bouncing boundary conditions increments are added to h and the interface grows, 
see also Ref. @. 

From an analysis of Eqs. (|4.1[) and (|4.2[) it follows that the growth modes with positive charge exerts an attraction 
of the other modes, whereas the negatively charged growth modes repel the other modes. Although our analysis so far 
only applies to a dilute network, it follows tentatively from a numerical simulations of Eqs. (|4. 1[) and (|4.2p that the 
modes form dipoles and that the long time stable network configuration is composed of a gas of propagating dipoles, 
i.e. the pairing of monopole growth modes with opposite charges. 

In terms of the slope field a dipole mode with charges k and —k has the form Udi P ~ (fc/fco)[(r — vt — r\)/\f — vt — 
f\\ — [f — vt — r2)/\r — vt — T2^, v = A(fc/fco)(r*2 — ^i)/|^2 — The mode propagate with velocity v and carries energy, 
momentum and action, E^i-p oc {v 2 /k 2 )k 4 ~ d , Hdi P oc (v/k 2 )k 3 ~ d , and Sdip oc (y 2 jk^)k 4 ~ d T . In terms of the height 
field hdi p ~ (k/ko)[\r — vt — r*i| — \r— vt — rb |] the mode in depicted in Fig. [4] Asymptotically the height field is flat 
corresponding to a vanishing slope field. The dipole mode corresponds to a propagating local defect or deformation. 
In Id the dipole mode corresponds in the slope field to a matched right hand and left hand domain wall propagating 
across the system, see Fig. S) 



The scaling issues for the KPZ equation remain unsettled except in Id where the fluctuation-dissipation theorem 
(|1.7p and the scaling law (|1.6p give access to the strong coupling fixed point with scaling exponents f = 1/2 and 
= 3/2, see Fig. [TJ In higher d there has been many attempts to access the strong coupling fixed point both on the 
basis of DRG, mode coupling, directed polymers (DP), and numerically; however, the strong coupling features remain 
elusive Q. 

The present weak noise method is not a scaling approach but rather a many body description of a growing interface. 
Nevertheless, the method allows a discussion of some of the scaling features. Since the scheme is consistently Galilean 
invariant the scaling law (|1.6p is automatically obeyed. The roughness exponent £ is associated with the static 
correlations (hh)(r) = J Y[dhh(r)h(0)Po(h) oc and requires the static distribution Po(h), Po(h) oc exp[— So(h)/A], 
So(h) = limr_ +00 S(hi — > h,T). The stationary distribution is associated with the zero-energy stationary manifold 
which in general is difficult to identify for a system with many degrees of freedom. 

In Id in terms of the slope field u = Vh and associated noise field p the Hamiltonian takes the form H = 
J dx p[vV 2 u + AmVu — (1/2) V 2 p] and it follows that the stationary zero-energy manifold is given by p — 2vu, yielding 
a Hamiltonian density as a total differential. Correspondingly, the stationary action So oc j dxdtpdu / dt = v J dx(Vh) 2 
and the distribution Pq(u) oc exp[— So/A] in accordance with (|1.7p . implying £ = 1/2. Since Galilean invariance is 
built in the exponent z = 3/2 follows automatically. However, we can also infer z = 3/2 from an independent 
argument based on the dispersion law for the low lying gapless excitations. The elementary excitations are the right 
hand and left hand domain walls. A composite quasi particle or dipole mode satisfying the boundary condition of 
vanishing slope can be constructed by paring two domain walls. The dipole mode propagates with energy E oc u 3 
and momentum n oc u 2 , where u is the dipole amplitude. Using the analogy between the stochastic formulation and 
quantum mechanics, i.e., the canonical quantization of the weak noise scheme with A as an effective Planck constant, 



V. 



SCALING AND UPPER CRITICAL DIMENSION 
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and using the spectral representation (uu)(x,t) = J dHF(U) exp(— Et/A + iHx/A), F(H) is a form factor, we note 
that the bottom of the gapless dipole energy spectrum E oc II 3 / 2 implies (uu)(x, t) oc G(t/x z ), where z = 3/2. 

In higher d we have not been able to identify the stationary zero-energy manifold and the exponent £. Also, in order 
to determine z, for example from the form of the time dependent correlations (hh){f, £), we need both the transition 
probabilities and the stationary distribution, i.e., by definition (hh)(f,t) = J Yl dhidh,2P{h\{0) — ► h2(r),t)Po(hi). In 
the weak noise scheme P{h\ — > h 2 ,t) oc exp[— S(hi — > h,2,t)/A] which requires a detailed analysis of the dynamical 
network and the associated action. 

In the dipole sector we can, however, present some preliminary scaling results. Since the propagating dipole mode 
according to (|3 .6[) carries action 5* oc k A ~ d T and propagates with velocity v oc k the mode moves the distance L = vT in 
time T. Expressing S in the form S oc L 4 ~ d /T 3 ~ d we infer the single dipole distribution P oc exp[— const. L 4 ~ d /T 3 ~ d ] 
and, correspondingly, the dipole mean square displacement (5L 2 ) oc T 2H with Hurst exponent H = (3 — d)/(4 — d); 
note that the dynamic exponent z = H^ 1 = (4 — d)/(3 — d). In the stochastic representation the dipole mode thus 
performs anomalous diffusion. In d = we have H = 3/4, z = 4/3, in agreement with a formal DP result [U, @]. 
In d — 1 we obtain H = 2/3 and the exact result z — 3/2, i.e., the dipole modes exhaust the spectrum. In d = 2 
we have if = 1/2 and z = 2, i.e., ordinary diffusion. In = 3 we obtain H ~ and z = oo, the dipole mean 
square displacement falls off logarithmically (SL 2 ) oc InT, however, z diverges at variance with accepted DRG and 
DP results. In d = 4 we have H = — oo and the mean square displacement is arrested. Below d = 2 (the lower critical 
dimension) the dipole modes superdiffuse, above d = 2 we have subdiffusion. These scaling results only refer to the 
dipole sector. 

The last issue is the much discussed upper critical dimension for the KPZ equation Q. The weak noise approach 
allows a non scaling argument for the existence of a critical dimension. Above d = 4 the negative growth mode as 
a bound state solution to the non linear Schrodinger equation (|3.4p ceases to exist. This implies that the dynamical 
network representation of a growing interface ceases to be valid. This result follows from a numerical analysis of the 
non linear Schrodinger equations (13. 4| but can also be inferred by an algebraic proof based on Derrick's theorem. First, 
introducing K = (1/2) J d d x(Vw) 2 , N = J d d xw 2 , and I = J d d xw 4 we infer from (|3~4"1) the identity -2K = k 2 N-k^I. 
Second, deducing (|3.4[) from a variational principle SF/Sw = with F = K + (\/2)k 2 N — (fc 2 /4)I and performing a 
constrained minimization w(f) — > w(/j,r), K — > /.i d ~ 2 K, N — > /jrN, I — > /j. d I and 5F/5/i\^ = i = we have the second 
identity (d - 2)K + (k 2 /2)dN - (k 2 /4)dl = 0. Finally, requiring N,I > the identities imply d < 4, q.e.d. 

VI. SUMMARY AND CONCLUSION 

In this paper we have presented a short review of a recently developed asymptotic weak noise approach to the 
Kardar-Parisi-Zhang equation. The scheme provides a many body description of a growing interface in terms of a 
dynamical network of growth modes. The growth modes are the elementary building blocks and their propagation 
accounts for the kinetic growth. Kinetic transitions are determined by an associated dynamical action, replacing the 
customary free energy landscape. Superimposed on the network is a gas of diffusive modes. In Id the dispersion laws 
delimit the universality classes: In the KPZ case the gapless domain wall modes yield z = 3/2, the diffusive modes 
being subdominant; in the EW case the domain walls are absent and the gapless diffusive modes yield z = 2. In 
higher d the scaling results based on the weak noise method are still subject to scrutiny. Finally, we mention that 
the weak noise method has also been applied to the noise-driven Ginzburg-Landau equation, a finitc-time-singularity 
model, and DNA bubble dynamics @. 
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